THE WAVE MECHANICS OF LARGE-SCALE STRUCTURE 



PETER COLES 

School of Physics & Astronomy, 
University of Nottingham, 
University Park, 
Nottingham NG7 2RD, 
United Kingdom 



Abstract. I review the basic "gravitational instability" model for the 
growth of structure in the expanding Universe. This model requires the 
existence of small initial irregularities in the density of a largely uniform 
universe. These grow through linear and non-linear stages to form a complex 
network of clusters, filaments and voids. The dynamical equations describ- 
ing the evolution of a self-gravitating fluid can be rewritten in the form 
of a Schrodinger equation coupled to a Poisson equation determining the 
gravitational potential. This approach has a number of interesting features, 
many of which were pointed out in a seminal paper by Widrow & Kaiser 
(1993). I argue that this approach has the potential to yield useful analytic 
insights into the dynamical growth of large-scale structure. As a particular 
example, I show that this approach yields an elegant reformulation of an 
idea due to Jones (1999) concerning the origin of lognormal intermittency 
in the galaxy distribution. 



1. Introduction 

The local Universe displays a rich hierarchical pattern of galaxy clustering 
that encompasses a large range of length scales, culminating in rich clus- 
ters and superclusters. The early Universe, however, was almost smooth, 
with only slight ripples seen in the cosmic microwave background radiation. 
Models of the evolution of structure link these observations through the ef- 
fect of gravity, because the small initially overdense fluctuations accrete 
additional matter as the Universe expands. During the early stages, the 
ripples evolve independently, like linear waves on the surface of deep water. 
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but as the structures grow in mass, they interact with other in non-hnear 
ways, more Uke nonhnear waves breaking in shallow water. 

The linear theory of perturbation growth is well-established, but the 
non-linear regime is much more complicated and generally not amenable to 
analytic solution. Numerical A^-body simulations have led the way towards 
an understanding of strongly developed clustering, but simulating a thing is 
not quite equivalent to understanding it. In this lecture, therefore, I sketch 
out some of the analytic methods that can be used to study non-linear 
clustering. I focus in particular on a novel approach based on the description 
of density fluctuations using quantum mechanics; see also Coles (2002). 



2. Cosmological Structure Formation 

The Big Bang theory is built upon the Cosmological Principle, a symmetry 
principle that requires the Universe on large scales to be both homoge- 
neous and isotropic. Space-times consistent with this requirement can be 
described by the Robertson-Walker metric 

where k is the spatial curvature, scaled so as to take the values or ±1. 
The case k = represents flat space sections, and the other two cases are 
space sections of constant positive or negative curvature, respectively. The 
time coordinate t is called cosmological proper time and it is singled out 
as a preferred time coordinate by the property of spatial homogeneity. The 
quantity a{t), the cosmic scale factor, describes the overall expansion of the 
universe as a function of time. If light emitted at time t^ is received by an 
observer at to then the redshift z of the source is given by 

a{te) 

The dynamics of an FRW universe are determined by the Einstein gravi- 
tational field equations which become 

- = SttGp ^ + A, 3 

a J a"' 

a 47rG / p\ A 

These equations determine the time evolution of the cosmic scale factor 
a{t) (the dots denote derivatives with respect to cosmological proper time 
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t) and therefore describe the global expansion or contraction of the universe. 
The behaviour of these models can further be parametrised in terms of the 
Hubble parameter H = (a/ a) and the density parameter ft = SirGp/SH"^ , 
a suffix representing the value of these quantities at the present epoch 
when t = to- 

In order to understand how of large-scale structure arises, it is best to 
begin with the standard fluid-based approach to structure growth. In the 
standard treatment of the Jeans Instability one begins with the dynamical 
equations governing the behaviour of a self-gravitating fluid. These are the 
Euler equation 

4t^ + (v-V)v + -Vp + V</. = 0; (6) 
ot p 

the continuity equation 

% + +V(pv) = , (7) 

expressing the conservation of matter; and the Poisson equation 

V^(/) = AttGp , (8) 

describing Newtonian gravity. If the length scale of the perturbations is 
smaller than the effective cosmological horizon dn = c/H, a Newtonian 
treatment of cosmic structure formation is still expected to be valid in 
expanding world models. In an expanding cosmological background, the 
Newtonian equations governing the motion of gravitating particles can be 
written in terms of 

X = r/a{t) (9) 

(the comoving spatial coordinate, which is fixed for observers moving with 
the Hubble expansion), 

v = r-Hr = ax (10) 

(the peculiar velocity field, representing departures of the matter motion 
from pure Hubble expansion), p(x, t) (the matter density), and ^(x, t) (the 
peculiar Newtonian gravitational potential, i.e. the fluctuations in potential 
with respect to the homogeneous background) determined by the Poisson 
equation in the form 

Vx^^ = AnGa^ip - po) = AnGa^pod. (11) 

In this equation and the following the suffix on indicates derivatives with 
respect to the new comoving coordinates. Here po is the mean background 
density, and 

6 ^ (12) 
po 
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is the density contrast. Using these variables the Euler equation becomes 

^igl + (v . Vx)v = -^VxP - Vx0 . (13) 

The first term on the right-hand-side of equation (13) arises from pressure 
gradients, and is neglected in dust-dominated cosmologies. Pressure effects 

may nevertheless be important in the the (collisional) baryonic component 
of the mass distribution when nonlinear structures eventually form. The 
second term on the right-hand side of equation (13) is the peculiar gravi- 
tational force, which can be written in terms of g = — Vx</'/a, the peculiar 
gravitational acceleration of the fluid clement. If the velocity flow is irrota- 
tional, V can be rewritten in terms of a velocity potential (py'. 

V = -Vxc/'^/a. (14) 

This is expected to be the case in the cosmological setting because (a) 
there are no sources of vorticity in these equations and (b) vortical pertur- 
bation modes decay with the expansion. We also have a revised form of the 
continuity equation: 

^ + 3Hp + - Vx(pv) = . (15) 
ot a 

In order to understand how structures form we need to consider the 
difficult problem of dealing with the evolution of inhomogeneities in the 
expanding Universe. We are helped in this task by the fact that we expect 
such inhomogeneities to be of very small amplitude early on so we can adopt 
a kind of perturbative approach, at least for the early stages of the problem. 
The procedure is to linearise the Euler, continuity and Poisson equations by 
perturbing physical quantities defined as functions of Eulcrian coordinates, 
i.e. relative to an unperturbed coordinate system. Expanding p, v and cp 
perturbatively and keeping only the first-order terms in equations (13) and 
(15) gives the linearised continuity equation: 

| = -iv..v, (16) 

ot a 

which can be inverted, with a suitable choice of boundary conditions, to 
yield 

The function / ~ ^^o^; this is simply a fitting formula to the full solution 
(Peebles 1980). The linearised Euler and Poisson equations are 



WAVE MECHANICS OF LARGE-SCALE STRUCTURE 



5 



(19) 



|v|, 1^1, \6\ <C 1 in equations (17), (18) & (19). From these equations, and if 
one ignores pressure forces, it is easy to obtain an equation for the evolution 
of S: 



For a spatially flat universe dominated by pressureless matter, po(*) = 
l/GirGt^ and equation (20) admits two linearly independent power law 
solutions d{x,t) = D±{t)5{x), where D^{t) oc a{t) oc i^/^ is the growing 
mode and -D-(i) oc is the decaying mode. 

The above considerations apply to the evolution of a single Fourier mode 
of the density field (5(x, t) = D^{t)S{x.). What is more likely to be relevant, 
however, is the case of a superposition of waves, resulting from some kind 
of stochastic process in which he density field consists of a superposition 
of such modes with different amplitudes. A statistical description of the 
initial perturbations is therefore required, and any comparison between 
theory and observations will also have to be statistical. Many versions of 
the inflationary scenario for the very early universe (Guth 1981; Guth & Pi 
1982; Brandenberger 1985) predict the initial density fluctuatiosn to take 
the form of a Gaussian random field in which the initial Fourier modes of 
the perturbation fleld have random phases. 

3. Nonlinear Gravitational Instability 

The linearised equations of motion provide an excellent description of grav- 
itational instability at very early times when density fluctuations are still 
small (5 -C 1). The linear regime of gravitational instability breaks down 
when S becomes comparable to unity, marking the commencement of the 
quasi-linear (or weakly non-linear) regime. During this regime the density 
contrast may remain small (6 < 1), but the phases of the Fourier compo- 
nents (5k become substantially different from their initial values resulting 
in the gradual development of a non-Gaussian distribution function if the 
primordial density fleld was Gaussian. In this regime the shape of the power- 
spectrum changes by virtue of a complicated cross-talk between diflFerent 
wave-modes. Analytic methods are available for this kind of problem (Sahni 
&; Coles 1995), but the usual approach is to use A'^-body experiments for 
strongly non- linear analyses (Davis et al. 1985; Jenkins et al. 1998). 

Further into the non-linear regime, bound structures form. The baryonic 
content of these objects may then become important dynamically: hydro- 
dynamical effects (e.g. shocks), star formation and heating and cooling of 
gas all come into play. The spatial distribution of galaxies may therefore 
be very different from the distribution of the (dark) matter, even on large 
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scales. Attempts are only just being made to model some of these processes 
with cosmological hydrodynamics codes, such as those based on Smoothed 
Particle Hydrodynamics (SPH; Monaghan 1992), but it is some measure of 
the difficulty of understanding the formation of galaxies and clusters that 
most studies have only just begun to attempt to include modelling the de- 
tailed physics of galaxy formation. In the front rank of theoretical efforts 
in this area are the so-called semi-analytical models which encode simple 
rules for the formation of stars within a framework of merger trees that 
allows the hierarchical nature of gravitational instability to be explicitly 
taken into account. 

Perturbation theory fails when nonlinearities develop but it is impor- 
tant to stress that the fluid treatment is intrinsically approximate anyway. 
A fuller treatment of the problem requires a solution of the Boltzmann equa- 
tion for the full phase-space distribution of the system /(x, v, t) coupled 
to the Poisson equation (8) that determines the gravitational potential. In 
cases where the matter component is collisionless, the Boltzmann equation 
takes the form of a Vlasov equation: 



The fluid approach outline above can only describe cold material where the 
velocity dispersion of particles is negligible. But even if the dark matter is 
cold, there may be hot components of baryonic material whose behaviour 
needs also to be understood. Moreover, the fluid approach assumes the 
existence of a single fluid velocity at every spatial position. It therefore 
fails when orbits cross and multi-streaming generates a range of particle 
velocities through a given point. 

Fortunately the formation of these structural elements can also be un- 
derstood using simple models, especially that of Zel'dovich (1970). This 
approximation actually predicts that the density in certain regions - called 
caustics - should become infinite, but the gravitational acceleration caused 
by these regions remains finite. Of course, in any case one cannot jus- 
tify ignoring pressure when the density becomes very high, for much the 
same reason as we discussed above in the context of spherical collapse: one 
forms shock waves which compress infalling material. At a certain point 
the process of accretion onto the caustic will stop: the condensed matter 
is contained by gravity within the final structure, while the matter which 
has not passed through the shock wave is held up by pressure. It has been 
calculated that about half the material inside the original fluctuation is 
reheated and compressed by the shock wave. An important property of the 
structures which thus form is that they are strongly unstable to fragmen- 
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tation. In principle, therefore, one can generate structure on smaller scales 
than the pancake. 

I will now describe the Zel'dovich approximation in more detail, and 
show how it can follow the evolution of perturbations until the formation of 
pancakes. Imagine that we begin with a set of particles which are uniformly 
distributed in space. Let the initial (i.e. Lagrangian) coordinate of a particle 
in this unperturbed distribution be q. Now each particle is subjected to a 
displacement corresponding to a density perturbation. In the Zel'dovich 
approximation the Eulerian coordinate of the particle at time t is 

r(t,q) = a(i)[q-6(i)Vq$o(q)], (22) 

where r = a(t)x, with x a comoving coordinate, and we have made a{t) 
dimensionless by dividing throughout by a{ti), where ti is some reference 
time which we take to be the initial time. The derivative on the right hand 
side is taken with respect to the Lagrangian coordinates. The dimensionless 
function b{t) describes the evolution of a perturbation in the linear regime, 
with the condition b{ti) = 0, and therefore solves the equation 

b + 2(^^^b-4TrGpb = ■ (23) 

cf. equation (20). For a flat matter-dominated universe we have b oc t^/^ 
as before. The quantity $o(q) is proportional to a velocity potential, of 
the type introduced above, i.e. a quantity of which the velocity field is the 
gradient: 

dr dx 
Hr = a — 

dt dt 

this means that the velocity field is irrotational. The quantity ^*o(ct) is 
related to the density perturbation in the linear regime by the relation 



V = — -Hr = a— = -a6Vq$o(q); (24) 



5 = bVl^o, 

which is a simple consequence of Poisson's equation. 

The Zel'dovich approximation is a linear approximation with respect to 
the particle displacements rather than the density, as was the linear solution 
we derived above. It is conventional to describe the Zel'dovich approxima- 
tion as a kind of first-order Lagrangian perturbation theory, while what 
we have dealt with so far for 5{t) is a first order Eulerian theory. We have 
also assumed that the position and time dependence of the displacement 
between initial and final positions can be separated. Notice that particles in 
the Zel'dovich approximation execute a kind of inertial motion on straight 
line trajectories. 
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The Zel'dovich approximation, though simple, has a number of interest- 
ing properties. First, it is exact for the case of one dimensional perturbations 
up to the moment of shell crossing. As wc have mentioned above, it also 
incorporates irrotational motion, which is required to be the case if it is 
generated only by the action of gravity (due to the Kelvin circulation theo- 
rem). For small displacements between r and a(i)q, one recovers the usual 
(Eulerian) linear regime: in fact, equation (22) defines a unique mapping 
between the coordinates q and r (as long as trajectories do not cross); this 
means that p{r,t)d^r = {p{ti))d^q or 

where |J(r, t)| is the determinant of the Jacobian of the mapping between 
q and r: dr/dq. Since the flow is irrotational the matrix J is symmetric 
and can therefore be locally diagonalised. Hence 

p{r,t) = {p{t))l[[l + b{t)a,{ci)]-': (26) 

i=l 

the quantities 1 + b{t)ai are the eigenvalues of the matrix J (the ai are 
the eigenvalues of the deformation tensor). For times close to ti, when 
\b{t) ai\ <^ 1, equation (26) yields 

6 -{ai + a2 + a3)b{t), (27) 

which is just the law of perturbation growth in the linear regime written a 
different way. 

At some time tgc, when b{tsc) = —^/oij-, an event called shell-crossing 
occurs such that a singularity appears and the density becomes formally 
infinite in a region where at least one of the eigenvalues (in this case aj) is 
negative. This condition corresponds to the situation where two points with 
different Lagrangian coordinates end up at the same Eulerian coordinate. 
In other words, particle trajectories have crossed and the mapping between 
Lagrangian and Eulerian space is no longer unique. A region where the 
shell-crossing occurs is called a caustic. For a fluid element to be collapsing, 
at least one of the aj must be negative. If more than one is negative, then 
collapse will occur first along the axis corresponding to the most negative 
eigenvalue. If there is no special symmetry, one therefore expects collapse to 
be generically one-dimensional, from three dimensions to two. Only if two 
(or three) negative eigenvalues, very improbably, are equal in magnitude 
can the collapse occur to a filament (or point). One therefore expects the 
formation of flattened structures to be the generic result of such collapse. 
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This is in accord with the classic work of Lin, Mestel k. Shu (1965) who 
showed that, for a generic triaxial perturbation, the collapse is expected to 
occur not to a point, but to a flattened structure of quasi-two dimensional 
nature. The usual descriptive term for such features is pancakes. 

The Zel'dovich approximation matches very well the evolution of den- 
sity perturbations in full iV-body calculations until the point where shell 
crossing occurs (Coles, Melott &: Shandarin 1993). After this, the approx- 
imation breaks down completely. Particles continue to move through the 
caustic in the same direction as they did before. Particles entering a pan- 
cake from either side merely sail through it and pass out the opposite side. 
The pancake therefore appears only instantaneously and is rapidly smeared 
out. In reality, the matter in the caustic would feel the strong gravity there 
and be pulled back towards it before it could escape through the other side. 
Since the Zel'dovich approximation is only kinematic it does not account 
for these close-range forces and the behaviour in the strongly non-linear 
regime is therefore described very poorly. Furthermore, this approximation 
cannot describe the formation of shocks and phenomena associated with 
pressure. 

Attempts to understand properties of non-linear structure using the 
fluid model therefore resort to further approximations (Sahni & Coles 1995) 
to extend the approach beyond shell-crossing. One relatively straightfor- 
ward way to extend the Zel'dovich approximation is through the so-called 
adhesion model (Gurbatov, Saichev & Shandarin 1989). In this model one 
assumes that the particles stick to each other when they enter a caustic 
region because of an artificial viscosity which is intended to simulate the 
action of strong gravitational effects inside the overdensity forming there. 
This 'sticking' results in a cancellation of the component of the velocity of 
the particle perpendicular to the caustic. If the caustic is two-dimensional, 
the particles will move in its plane until they reach a one-dimensional 
interface between two such planes. This would then form a filament. Mo- 
tion perpendicular to the filament would be cancelled, and the particles 
will flow along it until a point where two or more filaments intersect, thus 
forming a node. The smaller is the viscosity term, the thinner will be the 
sheets and filaments, and the more point-like will be the nodes. Outside 
these structures, the Zel'dovich approximation is still valid to high accuracy. 
Comparing simulations made within this approximation with full iV-body 
calculations shows that it is quite accurate for overdensities up to (5 ~ 10. 

The spatial distribution of particles obtained using the adhesion ap- 
proximation represents a sort of "skeleton" of the real structure: non-linear 
evolution generically leads to the formation of a quasi-cellular structure, 
which is a kind of "tessellation" of irregular polyhedra having pancakes for 
faces, filaments for edges and nodes at the vertices . This skeleton, however. 
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evolves continuously as structures merge and disrupt each other through 
tidal forces; gradually, as evolution proceeds, the characteristic scale of the 
structures increases. In order to interpret the observations wc have already 
described, one can think of the giant "voids" as being the regions internal 
to the cells, while the cell nodes correspond to giant clusters of galaxies. 
While analytical methods, such as the adhesion model, are useful for map- 
ping out the skeleton of structure formed during the non-linear phase, they 
are not adequate for describing the highly non-linear evolution within the 
densest clusters and superclusters. In particular, the adhesion model cannot 
be used to treat the process of merging and fragmentation of pancakes and 
filaments due to their own (local) gravitational instabilities, which must be 
done using full numerical computations. 



4. An Alternative Approach 

A novel approach to the study of collisionless matter, with applications 
to structure formation, was suggested by Widrow & Kaiser (1993). It iin- 
volves re-writing of the fluid equations given in Section 2 in the form of 
a non-linear Schrodinger equation. The equivalence between this and the 
fluid approach has been known for some time; see Spiegel (1980) for histor- 
ical comments. Originally the interest was to find a fluid interpretation of 
quantum mechanical effects, but in this context we shall use it to describe 
an entirely classical system. 

To begin with, consider the continuity equation and Euler equation for 
a curl-free flow in which v = V^), in response to some general potential V. 
In this case the continuity equation can be written 

^ + V-{pVcP) = 0. (28) 

It is convenient to take the flrst integral of the Euler equation, giving 

m(V«= = -^. (29) 

which is usually known as the Bernoulli equation. The trick then is to make 
a transformation of the form 

7p = a eyi]i{i(j) / u) , (30) 

where p = ipip* = = c?\ the wavefunction i/j{-K,t) evidently complex. 
Notice that the dimensions of u are the same as 4>, namely [L'^T~^]. After 
some algebra it emerges that the two equations above can be written in one 
equation of the form 

^^^ = -YV2^ + FV; + P^, (31) 
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where 

To accommodate gravity we need to couple equation (31) to the Poisson 
equation in the form 

= 47rGVV'* , (33) 

taking V to be (j). 

This system looks very similar to a Schrodingcr equation, except for the 
extra "operator" P. The role of this term becomes clearer if one leaves it 
out of equation (31) and works backwards. The result is an extra term in 
the Bernouilli equation that resembles a pressure gradient. This is often 
called the "quantum pressure" that arises when one tries to understand a 
quantum system in terms of a classical fluid behaviour. Leaving it out to 
model a collisionless fluid can be justified only if a varies only slowly on 
the scales of interest. On the other hand one can model situations in which 
one wishes to model genuine effects of pressure by adjusting (or omitting) 
this term in the wave equation. Widrow Sz Kaiser (1993) advocated simply 
writing 

ih^ = -—V^^ + m<t>{yi)^, (34) 

i.e. simply ignoring the quantum pressure term. In this spirit, the constant 
h is taken to be an adjustable parameter that controls the spatial resolution 
A through a dc Broglic relation A = fi/mv. In terms of the parameter v 
used above, we have v = Ti/m, giving the correct dimensions for Planck's 
constant. Note that the wavefunction if) encodes the velocity part of phase 
space in its argument through the ansatz 



VXx) = Vp(x)exp[i9(x/7i)], (35) 

where V0(x) = p(x), the local 'momentum field'. This formalism thus 
yields an elegant description of both the density and velocity fields in a 
single function. 

The approach outlined above is relatively new to galaxy clustering stud- 
ies, and many details still need to be worked out. One source of complexity 
arises when one places the system in an expanding context. To see what 
happens, let us define a scaled density x = p/Po = (1 + <^) and take J7 = 1. 
The continuity equation then becomes 

^ + V • (xV<^) = , (36) 

where the velocity potential ^ is now such that u = dx/dt = and a 
is the scale factor. It is convenient to take the first integral of the Euler 
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equation, giving 

I + ^(V,)- + (37) 

where 9 = 2^/3a'^H^ and is the gravitational potential. After some more 
algebra the system again becomes a Schrodinger-like wave equation, but in 
a rather than in t and using ilP' = such that 

w^ = -'^V^i^ + Vi^ + Pi^ , (38) 
oa I 

with V = <^-\-Q and P as before. 

5. The Origin of Spatial Intermittency 

Many types of non-linear system display a time-evolution characterized by 
the word "intermittency". While linear Gaussian processes involve fluctu- 
ations that are symmetric about their mean value, non-linear processes 
typically have highly skewed distributions. In the context of time series, in- 
termittent processes often have long quiescent periods punctuated by bursts 
of intense activity. In the spatial domain, intermittent processes are ones 
in which isolated regions of high density are separated by large voids; see 
Shandarin & Zel'dovich (1989). 

One particular aspect of galaxy clustering that has received some at- 
tention over the years has been the property that the one-point probal)ility 
distribution of density fluctuations p(p) appears to have a roughly log- 
normal form, i.e. logp has a roughly normal distribution (Coles &: Jones 
1991). The lognormal is a prime distribution producing intermittency, and 
was discussed in a pioneering paper by Kolmogorov (1962). Although in a 
qualitative sense the application of the concept of intermittency to large- 
scale structure seems plausible, a quantitative description of how it arises 
is not easy to obtain. Drawing on ideas discussed by Zel'dovich et al. (1985, 
1987), Jones (1999) suggested an analytical model for the cosmological con- 
text. On a simple level, this is quite easy to understand. If one has a linear 
process such that the output Y is constructed by co-adding a large number 
of independent contributing processes Xj, 

N 

Y = Y.X, (39) 

1=1 

as AT —> CO then the central limit theorem guarantees that Y is Gaussian 
as long as the Aj have finite variance. If one takes instead a multiplicative 
process of the form 

iV 

y = n^i (40) 

i=l 
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with the Xi stiU independent then the same theorem suggests that logY 
should be normal ss N ^ oo. Lognormal distributions consequently arise 
naturally in random multiplicative processes, such as those involving frag- 
mentation or coagulation. 

Whatever the details of its origin, it is now established that this property 
has an interesting connection with the scaling properties of moments of the 
probability of the distribution. Taking a generic random variable X, such 
that the distribution of X within cells of side L is denoted p{X\L), then 
the g-th moment at a given value of L is said to display scaling if 



This means that different powers q of the density field vary as a different 
power of the coarse-graining scale L. The function //(g) is called the inter- 
mittency exponent, and it can be extracted from observations. Jones, Coles 
&; Martinez (1992) showed that observations suggest a roughly quadratic 
dependence of fi(q) upon q and that this is related to the underlying near- 
lognormal form of the density fluctuations. A set displaying the form (41) 
is usually termed a multifractal; see Paladin &; Vulpiani (1987) for general 
discussion. 

We know that the distribution of density fluctuations is not exactly 
lognormal. The intermittency exponent can be written in the form 



where the Dq are scaling dimensions {D2 for example is the correlation 
dimension). We know that Dq (x q for a multifractal model whereas per- 
turbative methods suggest a simpler form of scaling such that Dq = Dq is 
constant, typical of a monofractal. 

One aspect of this is that the hierarchy of correlation functions that 
describe a lognormal distribution display Kirkwood (1935) scaling, while it 
appears from numerical A?^-body simulations that cosmological fluctuations 
display a different hierarchical form. For a discussion of the relationship 
between lognormal and hierarchical scaling, see Coles & Prenk (1991). 

It is within the overall framework of the fluid model that Jones (1999) 
sought to understand the observed intermittency of the large-scale structure 
of the Universe. Using the velocity potential introduced above, he first 
introduces an effective Bernoulli equation for the flow: 




(41) 



= -{q - i)A 



(42) 



(43) 



dt 2a2 



where 4> is the actual gravitational potential. This equation neglects terms 
involving pressure gradients as mentioned above. To cope with shell-crossing 
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events, Jones (1999) introduces an artificial viscosity v by adding a term 
to the right-hand-side of this equation: 



3t 2a^ -v^+^V^^.. (44) 

The viscosity is introduced to prevent the particle flow from entering the 
multi-stream region by causing particles to stick together at shell-crossing. 
This is also used in an approach called the adhesion approximation (Gur- 
batov, Saichev & Shandarin 1989). Using the Hopf-Cole transformation 
= —2i'\ogip and defining a scaled gravitational potential via (j) = 2ue 
we can write the Bernoulli equation as 

^ = vV\ + e{^)^. (45) 

This is called the random heat equation, because of the existence of the 
spatially-fluctuating potential term e(x). The gravitational potential changes 
very slowly even during nonlinear evolution (Brainerd, Scherrer & Villum- 
sen 1993; Bagla k, Padmanabhan 1994), so Jones (1999) assumes that it 
can be taken as constant and to be Gaussian distributed. An approximate 
scaling solution to (45) can then be found using a path integral adapted 
from that normally used in quantum physics (Feynman &; Hibbs 1965) ; see 
below for more details. In this approximation, the function 99 has a lognor- 
mal distribution (Coles & Jones 1991). We refer the reader to Jones (1999) 
for details; see also Zel'dovich et al. (1985, 1987). 

This model is one of the few attempts that have been made to under- 
stand the non-linear behaviour of the matter distribution using analytic 
methods. Although not rigorous it surely captures the essential factors in- 
volved. It does, however, suffer from a number of shortcomings. First, the 
approach does not follow material beyond the shell-crossing stage. Second, 
the viscosity v that is needed does not have properties that are very real- 
istic physically: it can depend neither on the density p nor the position x. 
Moreover, in the final analysis Jones (1999) takes the limit ^ 0, so it can- 
cels out anyway. One is tempted to speculate that its introduction may be 
unnecessary. Third, the function (/?(x, t) that emerges from equation (45) is 
not the desired density /9(x, t) nor does it bear a straightforward relation to 
the density. Finally, it is not clear how the motion of a coUisional baryonic 
component can be modelled within this framework. 

This formulation provides a useful complementary approach to many 
techniques, including AT-body methods. It also provides a new light with 
which to study the Jones (1999) model. Widrow k. Kaiser (1993) show 
using theoretical arguments and numerical simulations that this system 
allows accurate numerical evolution of the system beyond shell-crossing, so 
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it does not have the ad hoc construction needed by the Jones (1999) model 
to remedy this. 

Second, no artificial viscosity is required. Equations (45) and (34) are 
of the same form, apart from minor subtleties like the use of complex time 
coordinates. The potential term is easily understood in (34), and the wave- 
function ip now has a straightforward relationship to p. The upshot of this 
is that if one adopts the approximation of constant gravitational potential 
one can use the same path integral approach as described by Jones (1999). 
In a nutshell, given some initial wavefunction ^(x',f') one can determine 
the wavefunction at a subsequent time V'(x, t) using 



where the function K{x, t; x', t') involves a sum over all paths F connecting 
the initial and final states with t > t': 



where V is an appropriate measure on the set of classical space-time tra- 
jectories. For a particle moving in a potential T^(x, t) = m(l){-K, t) the action 
S for a given path F is given by 



Note the presence of the Gaussian field in equation (48) and hence in the 

exponential of the integrand on the right-hand-side of equation (47). To get 
an approximate solution to this system we can follow the same reasoning as 
Zeldovich et al. (1985, 1987) and Jones (1999), ignoring time-varying terms, 
using the Gaussian properties and counting the dominant contributions 
to the path integral to deduce that the integral produces a solution of 
lognormal form. This part of the argument is identical to that advanced by 
Jones, except that the solution is for ijj rather than ip and since p is 
then one directly obtains a lognormal distribution for the desired density 



It should be stressed that, although the present approach clearly pro- 
vides a more elegant formulation of the problem, the deduction of lognor- 
mality remains approximate; the lognormal is not the exact solution to the 
system to either Jones' equation (45) or to equation (34). How accurately 
this approximate form applies is open to doubt and will have to be checked 
by full numerical solutions. Interestingly, however, it is known to apply 
quite accurately in quantum systems such as disordered mesoscopic elec- 
tron configurations (Janssen 1998). As mentioned above, the Schrodinger 




(46) 




(47) 




(48) 



p(x,t). 
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approach yields a wavefunction t/; which is directly related to the parti- 
cle density p via p = IV'P- Such quantum systems also display lognormal 
scaling for properties such as the conductance, which depends on IV'I^- It 
is a property of the lognormal distribution that if a random variable X 
is lognormal, then so is X". In such systems the role of the gravitational 
potential (f) is played by a potential that describes the disorder of a solid, 
perhaps caused by the presence of defects. Such systems display localisation 
at low temperature which is similar in some ways to the original idea of 
Anderson location (Anderson 1958). The formation of strongly non- linear 
structures by gravity is thus directly analogous to the generation of localised 
wavefunctions in condensed matter systems. 

Finally, and perhaps most promisingly for future work, the equation 
(34) offers a relatively straightforward way of modelling the behaviour of 
collisional material. The addition to the potential of a term of the form 
(with K an appropriately-chosen constant), converts the original equation 
(34) into a nonlinear Schrodinger equation: 

Bib Ti^ 

= - — V^^ + m</.(x)V' + kIV'I V (49) 

(Sulem & Sulem 1999). This equation is now equivalent to those that de- 
scribe the flow of a barotropic fluid; see Spiegel (1980). This system can 
therefore be used to model pressure effects, which are otherwise only han- 
dled effectively using numerical methods such as smoothed-particle hydro- 
dynamical approximations (e.g. Monaghan 1992). In the context of quan- 
tum systems, the nonlinear term is used to describe the formation of Bose- 
Einstein condensates (e.g. Choi &; Niu 1999). 

6. Discussion 

In this lecture I have sketched out an approach to the study evolving cos- 
mological density fluctuations that relies on a transformation of the Vlasov- 
Poisson system into a Schrodinger-Poisson system. The transformation is 
not a new idea, but despite the efforts of Widrow & Kaiser (1993) it does 
not seem to be well known in the astronomical community. The immediate 
advantage of this new formalism is that it yields a rather more convincing 
approach to understanding the origin of spatial intcrmittcncy and approx- 
imate lognormality in the galaxy distribution than that offered by Jones 
(1999). It also makes a connection in the underlying physics with other 
systems that display similar phenomena. 

On the other hand, one must be aware of the approximations also inher- 
ent in the present approach. The Schrodinger equation is not exact, and its 
usefulness as an approximate tool is restricted by a number of conditions 
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outlined by Widrow & Kaiser (1993); see also Spiegel (1980). Furthermore, 
the lognormal solution of the system is a further approximation and may 
not be valid especially in the strongly- fluctuating limit. Although it neatly 
bypasses some of the problems inherent in the Jones (1999) model, the 
nonlinear wave equation is by no means easy to solve in general situations. 
Numerical methods will still have to be employed to understanding other 
aspects of the evolution of cosmic structure within this framework as indeed 
they are in other branches of physics. 

One particular issue worth exploring using this approach is to under- 
stand the limits of the approach in strongly non-linear situations. As it 
stands, the justification for the lognormal approximation arises from the 
weakly non-linear behaviour of collisionless matter moving in an almost 
constant potential field. Taking into account the expansion of the Universe, 
the changing gravitational potential, and the possible effects of matter pres- 
sure within in the action formalism may well reveal that a different form 
of hierarchical scaling pertains in the strongly non-linear regime. 

It is perhaps worth mentioning some specific ideas of things that could 
be done using this approach and for which there seem to be clear benefits. 

— Perturbation Theory. Standard perturbation methods do not guar- 
antee a density field that is everywhere positive. Re-casting cosmolog- 
ical perturbation theory in terms of ip, constructed so that ip"^ = p can 
remedy this. 

— Gas Pressure. Analytic techniques for modelling the effects of gas 
pressure are scarce, even in relatively simple systems such as Lyman-a 
absorption cloud (Matarrese & Mohayee 2002). The quantum pressure 
term (36) or alternative terms such as in (56) allow flexibility to model 
gas behaviour at least at a phenomenological level. 

— Shell-crossing. Methods such as the Zel'dovich approximation break 
down at shell-crossing, as described in Section 2.6. Although the simple 
ansatz I have used in this lecture docs assume a unique velocity at every 
fluid location, it is possible to construct more complex representations 
that allow for multi-streaming (Widrow &: Kaiser 1993). Note also that 
no singularities occur in the wavefunction at any time. 

— Reconstruction. It is interesting to speculate that it might be pos- 
sible to use the unitary structure of quantum mechanics in order to 
turn back the clock on evolved structure in order to reconstruct parts 
of the cosmic initial conditions. 

In general many techniques exist for studying the wave mechanics of 
disordered systems, such as the renormalization group and path-integral 
methods, few of which are used by cosmologists working in this area. It 
is to be hoped that the introduction of some of these methods may allow 
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better physical insights into the behaviour of non-hnear structure formation 
than can be found using brute-force A?^-body techniques. 
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